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We calculate the yield of high energy neutrinos produced in astrophysical sources with negligible 
magnetic fields varying their interaction depth from nearly transparent to opaque. We take into 
account the scattering of secondaries on background photons as well as the direct production of 
neutrinos in decays of charm mesons. If multiple scattering of nucleons becomes important, the 
neutrino spectra from meson and muon decays are strongly modified with respect to transparent 
sources. Characteristic for neutrino sources containing photons as scattering targets is a strong 
energy-dependence of the ratio 7?" of and Ue fluxes at the sources, ranging from R° = (P^/fj), ~ 
below threshold to ~ 4 close to the energy where the decay length of charged pions and kaons 
equals their interaction length on target photons. Above this energy, the neutrino flux is strongly 
suppressed and depends mainly on charm production. 
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I. INTRODUCTION 

Experimental high energy neutrino physics has become 
one of the most active areas of astroparticlc physics, offer- 
ing among others the prospect of identifying the sources 
of ultra- high energy cosmic rays 0. High energy neu- 
trinos from astrophysical sources are the decay products 
of secondary mesons produced by scattered high energy 
protons on background protons or photons. The classic 
example arc the so-called cosmogenic or GZK neutrinos 
produced in scatterings of extragalactic ultra-high energy 
cosmic rays on cosmic microwave photons during propa- 
gation 01 . Two different kind of bounds on high energy 
neutrino fluxes exist: The cascade or EGRET limit uses 
bounds on the diffuse MeV-GeV photon background to 
limit the energy transferred to electromagnctically inter- 
acting particles that are produced unavoidably together 
with neutrinos Q. The cosmic ray upper bounds of, 
e.g., Refs. use the observed ultra-high energy cos- 

mic ray flux to limit possible neutrino fluxes. The latter 
limit assumes that all neutrino sources are transparent to 
hadronic interactions and thus at least neutrons can es- 
cape from the source region without interactions. Drop- 
ping the assumption of transparent sources and bidding 
the acceleration region by sufficient material absorbing 
ultra-high energy cosmic rays allows one to avoid the cos- 
mic ray limits [a, S B ■ One might therefore speculate 
that large neutrino fluxes at high and ultra-high energies 
might be produced in opaque sources, as they are needed, 
e.g., to perform neutrino absorption spectroscopy |lOl |. 

In this work, we calculate the flux of high energy neu- 
trinos produced as secondaries in astrophysical sources. 
We consider sources with interaction depth ranging from 
nearly transparent to opaque. In contrast to most earlier 
investigations |^ 0, 0| , we put emphasis on sources 
with such high densities that magnetic fields can be ne- 
glected and multiple scatterings are important. As a 
result, the neutrino spectra from meson and muon de- 
cays are strongly modified with respect to transparent 



sources. Since in most astrophysical environments the 
depth for is above threshold much larger than for pp 
interactions, we consider photons in the source as target 
material. In Sec. 2, we discuss the relevant production 
processes and present the neutrino yields from a single 
source. We discuss also briefly the consequences of the 
derived energy dependence of the neutrino yields on the 
expected high energy neutrino fluxes from astrophysical 
sources. In Sec. 3, we examine the neutrino flavor com- 
position at the source as well as the effect of neutrinos 
oscillations on the flavor composition at the detector and 
find significant deviations from the canonical flavor ra- 
tio expected from pion decay. Finally, we summarize our 
results in Sec. 4. 



II. NEUTRINO PRODUCTION PROCESSES 
AND YIELDS 

A. Simulation and particle interactions 

We idealize a hidden neutrino source as an homoge- 
neous slab filled with photons in which high energy pro- 
tons are injected. The probability Af that a particle trav- 
els in the slab the distance AL without scattering or de- 
cay is given by 

where I1/2 and knt are its decay and interaction length, 
respectively. Their relative value determines the fate of 
the particle, either it decays or scatters. 

In our Monte Carlo simulation, we track explicitly all 
secondaries (N, tt*, K^, K1 5) for which the interaction 
rate is non-negligible compared to their decay rate. We 
determine the multiplicity and energy spectra of light sec- 
ondary particles produced in scattering processes with a 
modified version of SOPHIA ■ For the case of hadrons 
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containing charm quarks, we employ HERWIG to de- 
termine the differential energy spectra of prompt charm 
neutrinos, while we estimate the total cross section for 
charm production from Ref . |l6l | . To fix the pion-photon 
and kaon-photon interactions, we use the following sim- 
ple recipe |T^ : We determine first the maximum of the 
resonant production cross section in the Breit-Wigncr ap- 
proximation of p and K* , respectively, and compare it to 
p + 7 ^ A ^ all. Then we rescale the proton-photon 
cross section by the ratios of the maxima to obtain the 
pion-photon and kaon-photon cross sections in the non- 
resonant region. Finally, the interaction length Zjnt is 
calculated from 

K^t - ^2 r de' a,,{e') r (2) 

where T is the Lorcntz factor of the hadron in the lab 
frame, £th is the threshold energy of the relevant reaction 
in the rest system of the hadron, and n{e) the energy 
distribution of photons with energy e in the lab frame. 

B. Phenomenological characterization of tlie 
sources 

In the following, we do not consider any particular 
model of neutrino sources, but characterize instead the 
sources in a phenomenological way. This allows us to per- 
form a systematic analysis of different sources according 
to their thickness, from nearly transparent ones where 
the accelerated protons hardly scatter once with the sur- 
rounding material to completely opaque sources where 
neither high energy cosmic rays nor high energy photons 
can escape. In our idealization of a neutrino source as an 
homogeneous slab filled with photons, each source is fully 
characterized by its length L and the photon distribution 
n{e). Potential sources of photon backgrounds are man- 
ifold, but we restrict ourselves in this work to the most 
important example, a background of thermal photons at 
temperature T. 

A more severe restriction of the current work is that 
we consider only sources with negligible magnetic fields. 
Thus we require that i) energy losses due to synchrotron 
losses are much smaller than due to interactions and that 
ii) deflections of charged particles in magnetic fields are 
small compared to the size of the source. This limita- 
tion allows us to stress more clearly the salient features 
of opaque sources. In a subsequent work, we shall con- 
sider the generic case of sources with arbitrary interaction 
depth and magnetic fields [T8l |. 

For the discussion of our numerical results it is use- 
ful to introduce following Ref. ^3 several dimcnsionlcss 
quantities. The first one is the "interaction depth" of 
nucleons defined analogously to the optical depth as the 
ratio T = L /li^t of the size L of the source to the inter- 
action length /int of nucleons. This ratio determines, if 
most of the protons leave the source without interactions 
(r ^ 1 or "transparent source" ) , or if multiple-scattering 



of nucleons is important and mesons are efficiently pro- 
duced (t ^ 1 or "opaque source"). For the illustration 
of our numerical results, we determine in the following 
T via /int = l/('^70') with a = 0.2 mb as reference cross 
section. 

As second parameter we introduce the dimensionless 
energy x = ELu/nip, where uj — 1.6T denotes the energy 
of the maximum of the Planck distribution and nip the 
proton mass. This parameter allows one to express the 
different yields in the form of a universal function, which 
in the case of transparent sources does not explicitly de- 
pend on the temperature. Finally, it is useful to introduce 
the critical dimensionless energy Xa = E^-cIjJ /rn^, where 
-Bcr = m/(cm£Q) is defined as the energy at which the de- 
cay length £i^2 = = T^o of a meson with life-time 
Ti/2 and mass m is equal to its interaction length ^int. 
Using charged pions as reference particles, x^ depends 
numerically as Xcr « 2x 10^°/(r/K)^ on the temperature 
of thermal photons. Thus, the parameter t determines 
the fraction of nucleons that scatter on photons, while the 
parameter x controls if mesons mainly decay (x <C Xcr) 
or scatter (a; ^ Xcr)- 

Before we discuss the details of the different spec- 
tra, we briefly comment on the main characteristics of 
hadron-photon interactions. In the upper panel of Fig. ^ 
we show the interaction length of protons, pions and 
kaons in a thermal bath of photons with temperature 
T = K together with their decay length as func- 
tion of the dimensionless energy x (bottom axis) or their 
energy (upper axis). The fate of a hadron is character- 
ized by two important energies: The threshold energy at 
a^th ~ 0.1, below which photo- meson reactions are expo- 
nentially suppressed, and the critical energy Xcr, above 
which most mesons will scatter before decaying. If the 
decay length and the interaction length cross inside the 
source, I1/2 = lint < L, the effect of multiple scatter- 
ing has to be taken into account. Main consequence is a 
strong suppression of the high energy neutrino flux pro- 
duced in meson decays for x > Xcr- The shorter life-time 
of kaons with respect to pions as well as their larger mass 
leads to a shift of their critical energy, i.e. a;^ > x^lr- 

In the lower panel of Fig. ^ we show the interaction 
and decay length of hadrons as a function of x for sources 
with different temperature T. As expected from the tem- 
perature dependence of ^int and the definition of x, the 
critical energy Xcr does not only depend on the meson 
considered but also on the temperature: Xcr scales sim- 
ply as Xci- oc l/T^ as long as Xcr is above threshold, 
Xcr ^ '--^ 0.1. Thus, multiple scattering starts to be- 
come important at lower energies as the source temper- 
ature increases. On the other hand, the dimensionless 
threshold energy xth is nearly independent of T. There- 
fore in the case of opaque sources the range with unsup- 
pressed neutrino fluxes, 0.1 ~ xth < x < 
shrinks for increasing T and becomes a narrow peak 
around xth for T > 5 x 10^ K. On the other hand, large 
neutrino fluxes at energies high enough to perform neu- 
trino spectroscopy, i?cr ^ 10^^, can be obtained only in 
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C. Neutrino yields and fluxes from a single source 
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FIG. 1: (Color online) Top: Interaction length (solid lines) 
for proton in black, charged pions in red (gray) and charged 
kaons in blue (dark gray) with thermal photons at a tempera- 
ture T = 10* K, together with the decay length (dashed lines) 
for TT*, and Kg. Bottom: Interaction length (solid lines) 
for proton and decay length (dashed lines) of charged pions 
for sources at T = 10'' K in black , T = lO'^ K in red (gray), 
and T — 10^ K in blue (dark gray). 



In order to analyze the effect of the thickness of the 
source on the neutrino spectrum we study the neutrino 
yield of a single source. The neutrino yield Y^{E), i.e. 
the ratio Y^{E) = (l)i,{E)/{T(j)p{E)) of the emitted neu- 
trino flux and the product of the depth r and the in- 
jected proton flux 0p, represents the number of neutrinos 
produced per injected proton with the same energy. For 
the energy spectrum of the injected protons we assume a 
power law dN/dE oc E~°' with a = 2.2, consistent with 
typical predictions for Fermi shock acceleration. For il- 
lustration wc choose the maximal energy of the initial 
proton spectrum as i?max — 10^"' eV. 

Let us first describe the different reactions contribut- 
ing to the neutrino flux. In Fig. |21 we show the var- 
ious contributions to for the case of an intermedi- 
ate value of the interaction depth, r = 0.4 j^^. The 
low energy tail of neutrinos arises from the decay of 
neutrons, n ^ p + + i>e, where the energy fraction 
transferred to neutrinos is on average only ~ 10""^. At 
all other energies, charged pions are produced most effi- 
ciently and their decay products provide the dominating 
contribution to the neutrino yield. The various decay 
channels of kaons contribute around 10% to the neutrino 
yields. The neutrino yield from prompt decays of charm 
is even smaller, because of the relative high mass of charm 
mesons resulting in a high energy threshold and small 
cross section for charm production. 




x=Ea)/ m„ 



sources at low temperatures, T < 5 x 10'^ K. If, addition- 
ally, one requires the source to be opaque so that the 
high-energy cosmic ray flux is suppressed, e.g. r > 10, 
then a second condition for a suitable source is a large 
extension, L > 3 x 10'^'^ / (T /K)^ cm. 



FIG. 2: (Color online) Neutrino yield Yi, as function of 2: = 
Euj/rrip from the decay of vr^ in magenta (light gray), 
in solid blue (dark gray), K'^ in dashed blue (dark gray), 
neutrons in green (very light gray), and charm mesons in red 
(gray) for a source with T = 10^ K and r — 0.4. 



Next wc discuss the effect of multiple scatterings on 
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the different neutrino yields. In the top panel of Fig. |31 
we show the yields from meson decays for a source with 



T 



10 K and three different sizes, L 
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(r = 0.04), 10" cm (t = 0.4) and lO^^ cm (r = 40), 
respectively. For r < 1, multiple scattering is negligible, 

oc r and the neutrino yields expressed as functions 
of X are thus independent from r. Moreover, the neu- 
trino yield from charm meson decays is clearly subdom- 
inant. Increasing the size of the source, multiple scat- 
tering cannot be neglected anymore and three different 
X regions can be distinguished. For x < x'^^, most pions 
decay before scattering and, as in the case of transparent 
sources, the neutrino yield is dominated by the contri- 
bution of charged pions. In the intermediate x range, 

^ ^ ^ ^CTi multiple pion scattering on photons be- 
comes effective and therefore neutrinos from kaon decays 
start to provide the most important contribution to the 
total neutrino yield |29l |. Going to even higher energies, 
x > , leads to a strong suppression of neutrino from 
these decays as both pions and kaons most often scatter 
before decaying. Hence at high energies, neutrinos from 
decays of charm mesons represent the main component of 
the total neutrino yield. Moreover, the maximum of the 
neutrino yield does not coincide with that of the trans- 
parent cases, because (jii, is not longer proportional to the 
interaction depth. 

In the bottom panel of Fig. |3 we show again neu- 
trino yields but for a source with T = 10^ K. If multiple 
scattering can be neglected, the neutrino yields are as ex- 
pected independent from T and r. For r > 1, we observe 
the same behavior of the neutrino yields from pion and 
kaon decays as in the upper panel, but the suppression 
of their fluxes starts already at lower energies. As antic- 
ipated in the discussion of Fig. 1, the range of energies 
where neutrino fluxes are unsuppressed strongly depends 
on the energy threshold Xth and the critical energy Xcr- 
So, whereas for a source at T = 10''' K the plateau visible 
in Fig. 121 for extends approximately over four orders 
of magnitude, in the case of a source at T = 10^ K it is 
reduced by two orders of magnitude. 

Finally, we want to illustrate how the high-energy sup- 
pression of neutrino fluxes depends on the interaction 
depth T, as well as on the density and the typical energy 
of photons. Figure 0] shows (unnormalized) fluxes for a 
source at T = 10^ K and different interaction depths, 
from 0.4 to 4 x 10"^. In the case of thin sources, r < a 
few, we observe that the final proton flux is only slightly 
distorted relative to the initial flux, and therefore pro- 
vides a non-negligible contribution to the observed UHE- 
CRs. The neutrino flux in these sources is basically pro- 
portional to the depth. For opaque sources the final 
proton flux becomes strongly suppressed above thresh- 
old, w niprriTr / {2e^) . For such high depths most 
pions and kaons scatter before decaying and as a con- 
sequence the flux of neutrinos becomes also suppressed 
at energies higher than Ecr- The only effect of further 
increasing t is an increase of the low energy neutrinos, 
as the high number of scatterings leads to more low en- 



E(eV) 

10'" 10'^ 10'^ 10^° 10^^ 

I iiiiii| I I iiiiii| I I iiiiii| I I iiiiii| I I iiiiii| I I iiiiii| I I a\ 
T=40 




E(eV) 



. 1 



10'^ 10'* 10'^ 10'' 10^' 



H 1 1 iiiiii| 1 1 iiiiii| 1 1 iiiiii| 1 1 iiiiii| 1 1 iiiiii| 1 1 iiiiii| 1 1 iiiiii| 1 1 iiiiii| 1 1 iiiiii| 1 1 m 

T=0,04 T=40 

T=0.4 




x=Ew/ rrip^ 



FIG. 3: (Color online) Neutrino yield Y^, as function of a: = 
Ecj/rrip (bottom axis) and E (upper axis) from pion magenta 
(light gray), kaon in blue (dark gray), and charm decays in red 
(gray) for r = 0.04, 0.4 and 40; upper panel for T = lO" K, 
lower panel for T = 10^ K. 



ergy mesons. Figure |21 shows (unnormalized) fluxes for 
an opaque source, r = 40, at T = 10^ K in the top and 
T = 10'' K in the bottom. As previously mentioned, the 
final proton flux is strongly suppressed above threshold, 
whereas the neutrino flux rises already at lower energies, 
because they carry away only a fraction of the proton en- 
ergy. The energy range where the neutrino flux is max- 
imal extends approximately from E^ to Ea, which for 
the source at T = 10^ K, goes roughly from 10^^ eV to 
10^" eV. This range is strongly reduced for increasing T 
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as Ecr (X 1/n oc l/T^, whereas i?th oc 1/T. Therefore, 
critical and threshold energies merge and only a a nar- 
row energy range around E 10^^ eV with unsuppressed 
neutrino flux remains for T = 10^ K, cf. the lower panel 
of Fig.m 
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FIG. 4: (Color online) Unnormalized fluxes of initial (dashed) 
and final (dotted) protons, as well as the total neutrino (solid), 
for a source at T = IC K and with interaction depths r = 0.4 
in magenta (light gray), 4 in blue (dark gray), 40 in red (gray) 
and 4 x 10^ in black. 



III. FLAVOR DEPENDENCE OF NEUTRINO 
YIELDS 

Both neutrino telescopes and extensive air shower ex- 
periments have some flavor discrimination possibilities. 
The case of neutrino telescopes is discussed for the ex- 
ample of ICECUBE in detail in Ref. . The long range 
of muons ensures that a muon track from charged- 
current interactions is always visible allowing the iden- 
tification of these events. By contrast, the charged- 
current interactions of j/g and Vr are — as long as the 
tau decay length is too short to be detectable — only 
distinguishable by the different muon content in elec- 
tromagnetic and hadronic showers. Therefore, these 
two flavors are in practice impossible to differentiate for 
E < 5 X 10^^ eV at a neutrino telescope. In the energy 
range 5 x 10"eV < < 2 x lO^^eV, the "double-bang" 
signature of gives some handle for the identiflca- 
tion of this flavor, while for higher energies at least one of 
the two shower is outside a 1 km^ detector. On the other 
hand, extensive air shower experiments have the poten- 
tial to identify double-bang events in the small energy 
interval 5 x 10"eV < < 2 x lO^^eV In summary. 
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FIG. 5: (Color online) Unnormalized fluxes of initial (dashed) 
and final (dotted) protons, as well as the total neutrino (solid), 
for a source with interaction depth r = 40 at T = 10* K (top) 
and T — 10^ K (bottom). The different contributions to the 
total neutrino fiux are also shown: pion in magenta (light 
gray), kaon in blue (dark gray), charm meson in red (gray), 
and neutron decay in green (very light gray). 



the main observable for the neutrino flavor composition 
is the ratio of track to shower events in a neutrino tele- 
scope, — (j)^/{(j)e + 4>t)i while only in a very small 
energy range all flavors can be distinguished. Addition- 
ally, extensive air shower experiments are sensitive to the 
fraction of tau events in all horizontal neutrino events, 
= 0^/(02 -|- (/)^), in a small energy window around 
10^8 eV. 



6 



In spite of these flavor discrimination possibilities, the 
potential of high energy neutrino observations for mixing 
parameter studies are only rarely discussed (for some re- 
cent examples see Ref . |2J, |2J| ) . One reason for this miss- 
ing interest is the small number of events expected in next 
generation experiments like ICECUBE or AUGER even 
in optimistic scenarios. Moreover, the neutrino flavor ra- 
tio from pion decay is 4>{ve) ■ ■ 4'{i^t) = 1:2:0 
before oscillations, resulting into (t>{ye) ■ '■ 4'{'^t) = 
1 : 1 : 1 at the detector quite independent from unknown 
details of the neutrino mixing matrix . Thus the max- 
imal mu-tau mixing together with the initial flavor ratio 
(pii^e) ■ ■ 4>{^t) = 1:2:0 seems to disfavor high en- 
ergy neutrino experiments as tools to probe only poorly 
known parameter as i?i3 or completely unknown ones as 
the octant of i?23 or the CP-violating phase 5c, v o f the 
neutrino mixing matrix. As discussed in Rcfs. |24L |25| , 
there exist however several examples of neutrino sources 
where at least in some energy range significant devia- 
tions from this canonical flavor ratio can be expected. 
Such a deviation does not only contain information about 
neutrino properties but also their sources. For instance, 
Ref. [2^ discussed the possibility to distinguish between 
pp or p7 sources of neutrinos measuring their flavor ra- 
tios. In the following, we will show that opaque sources 
arc characterized by an energy-dependent flavor ratio and 
thus the flavor ration encodes non-trivial information. 



A. Energy dependence of the neutrino flavor ratio 
at the source 

We start with an analysis of the expected neutrino fla- 
vor ratio W = Yy^jYv^ at the source. In Fig.|B| we show 
separately for each reaction contributing to the neu- 
trino yield in the case of a source with T = 10^ K. For 
definiteness we have set i?° = 0, if ^i,{x)l^p{E) < 10^^. 
The two panels show i?" for transparent source with in- 
teraction depth r = 0.04 (top) and r = 40 (bottom), 
respectively . 

In pion decays, all three neutrinos produced have very 
similar energy spectra and therefore is close to two. 
The most important decay channel of charged kaons is 
the same as the one of charged pions, fi + v^j^ 

2v^ + V(, + e. However, the mass difference between kaons 
and muons leads to very different energy distributions of 
the directly produced neutrinos. The neutrino spectra 
from kaons are dominated at low energies by neutrinos 
from muon decay, whereas at high energies z^^'s produced 
in the direct kaon decays dominate. As a consequence, 
there are three different ranges for the flavor ratio i?" of 
neutrinos produced by the charged kaon decays: At low 
energies the contribution of the direct is negligible. 
Thus EP approaches one at low energies, since the 
and Ve spectra from muon decay are similar. Near the 
threshold energy for kaon production, « 0.2, direct 
z/^'s start to dominate, leading to an increase of i?" to val- 
ues larger than two. The particular value of within 



this energy range depends on the initial kaon spectrum: 
the steeper it is, the larger the flavor ratio. As the energy 
increases the neutrino flavor ratio remains roughly con- 
stant until it gets closer to the high energy kinematical 
limit. At this point there are practically only z^^'s from 
the direct kaon decay and therefore BP increases fast. In 
the case of neutrinos produced in neutral kaon decays 
one expects a flavor ratio between two, corresponding to 
the contribution from (via 27r), and roughly 

one, due to the neutrino yield from the decay, mainly 
from Kl — *■ TT + e + Ue- Finally, concerning the neutri- 
nos generated in the decay of charm mesons one observes 
the same contribution from and t'^ , leading to a flavor 
ratio of one [s^ • 

The flavor ratio for an opaque source with r = 40 
(bottom panel) is only changed for neutrinos from kaon 
decays. At energies larger than x^' most of charged kaons 
scatter before decaying, leading to a steeper spectrum of 
the kaons which eventually decay. As a consequence at 
that point the value of i?° increases a second time before 
approaching the high energy limit. Roughly at the same 
energies the pions emitted by scatter before decaying. 
Therefore, in contrast to transparent sources, the flavor 
content of i^^'s at high energies is dominated by i/g and 
the ratio tends to values smaller than one. 

Let us now consider the total flavor ratio expected for 
different sources, shown in the Fig. [7| In the case of a 
transparent source the value of is small at low x be- 
cause almost all neutrinos are produced in the decay of 
neutrons, see Fig.|2| Once the contribution from charged 
pion decay becomes dominating, i?" approaches the stan- 
dard value two, as expected from the upper panel of 
Fig. |S1 In contrast to transparent sources, in the case 
of an opaque source, t = 40, R'^ is strongly energy- 
dependent. In particular, the flavor ratio reaches the 
value of 2 at lower x, because multiple scattering of nu- 
clcons leads to lower energies of the escaping neutrons. 
Moreover, there is a bump in the range 100 < x < 
and 1 < a: < 100 for the case T = 10'' K and 10^ K, 
respectively. The lower limit of theses ranges corre- 
sponds to the cross-over between pion and kaon domi- 
nance, while the upper limit corresponds to the cross- 
over between kaon and charm dominance. As can be in- 
ferred from the bottom panel of Fig.El the value of R'^ at 
the bump is determined by the interplay of the charged 
and neutral kaon neutrinos. Finally, in both cases the 
asymptotic value is R'^ ^ 1, which corresponds to the ra- 
tio expected when the decay of charm mesons dominate. 



B. Efftect of neutrino oscillations 

The neutrino spectra at the source discussed in the pre- 
ceding subsection are modulated by oscillations. There- 
fore the expected flavor ratios Ri at the Earth are differ- 
ent from the original ratios at the source. 

The neutrino fluxes arriving at the detector, (j>^ , can 
be written in terms of the initial fluxes 0q and the con- 



7 



E(eV) 



10^ 10 10" 10° 10^" 10 



iiij I iiiiiii| 



7T decay 
K"*" decay 
K° decay 
charm decay 




10"' 10"* 10'' 1 10' 10' 

x=Ea)/ trip' 

E(eV) 

10'' 10'" 10" 10'° 10'° 10" 



a: 



T 

9 

8 - 
7 
6 - 
5 
4 
3 - 



iiij I iiiiiii| 




7T decay 
K"*" decay 
K° decay 
charm decay 




10' 
x=Ea)/ m„ 



FIG. 6: (Color online) Flavor ratio R° of the neutrinos pro- 
duced in a source before they oscillate in terms of a; = Euj/rrip. 
Top: Specific flavor ratio for the neutrinos generated from the 
decay of tt^ in magenta (light gray), in solid blue (dark 
gray), Tv" in dashed blue (dark gray), and charm meson in 
red (gray) for a source with T = lO"* K and r = 0.04. In black 
dotted line is shown the total i?". Bottom: The same for a 
source with interaction depth t — 40. 



version probabilities Pa/s, 



if = ^ PaP<Pf3 = Paeffle + Pati4>ti. ■ 




(3) 



Since the interference terms sensitive to the mass split- 
tings Am^'s do not contribute, the conversion probabili- 
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FIG. 7: (Color online) Total flavor ratio for a source at T = 
10" K and T = 0.4 in dashed blue (dark gray), T = lO'' K 
and r = 40 in solid blue (dark gray), and T = 10'' K and 
r = 40 in solid red (gray). The main reactions giving rise to 
the neutrino signal are also shown. 



ties are simply 



P^p = <5„/3 - 2 ^ 5R(C/^^.C/^,.C/„,;7^,) , 
j>k 



(4) 



where U is the neutrino mixing matrix and Greek (Latin) 
letters are used as flavor (mass) indices. In Fig. |S1 we 
show the expected flavor ratio i?^ at the detector in the 
case of a source with T = 10^ K and r = 40 as a solid, 
red line for the best-fit point of the neutrino mixing 
rameters, sin^ §12 = 0.3, sin^ 1923 = 0.5 and -di^ = 
Various bands show the range of i?^ allowed if the mix- 
ing parameters are varied within their 95% C.L., while 
the unconstrained CP phase i5cp is varied in the whole 
possible range, Sep € [0 : tt]. Dominant uncertainty for 
the prediction of i?^ is the value of i923- 

We analyze next the dependence of the flavor ratio i?^ 
on the different neutrino mixing parameters in more de- 
tail. Expanding i?^ up to first order in ?9i3 and choosing 
for illustration i9i2 = n/6 gives 



R,.iR°, ^23, ^13, Sep) 



Pc 



RnP, 



Pee + RoPefi. + Per + RoPfiT 

- A(i?o, ^23) + 5(7?°, 7923) COS Scp^i3 + 0(1923). (5) 

The explicit expression for A and B are given in the 
Appendix. In Fig. O we show the dependence of the 
coefficients A (left panel) and B (right panel) on i?" and 
1923. While the parameter A is more sensitive to the value 
of if 7923 is in the second octant, sin^ t923 > 0.5, this 
behavior is opposite for B. 
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10"* 10"' 1 10' 10' 

x=Ea)/ trip' 

FIG. 8: (Color online) Flavor ratio i?^ at the Earth for a 
source with T = IQp K and r = 40. The solid red line (gray) 
corresponds to the best-fit point of the neutrino mixing pa- 
rameters, sin^ i?i2 = 0.3, sin^ i923 ~ 0.5 and dri = 0. In the 
colored areas mixing angles within 95% C.L. and a possible 
non-zero value of 5cp have been considered: iJia and 5c p in 
yellow (very light gray), plus i923 in blue (dark gray), and i9i2 
in green (light gray). 



Let us first analyze the dependence of on 'i?23. 
For the sake of simplicity, we assume i?i3 — 0, and 
thus i?^(i?°,?923,0,(Scp) = A[R^,d2z). The upper panel 
of Fig. ^1 shows the flavor ratio for different val- 
ues of sin^ i?23 . The possibility to extract information 
about the value of t?23 strongly depends on the energy 
range considered. At low energies, where the initial ra- 
tio is close to zero (neutron neutrinos), i?^ depends in- 
versely on sin^"i?23, varying between 0.15 and 0.35 for 
sin^ i?23 = 0.65 and 0.35, respectively. This large sen- 
sitivity of Rp, to the value of 7^23 makes the neutrinos 
emitted in neutron decay an alternative tool to study 
this mixing angle |25l |. At intermediate energies the sen- 
sitivity becomes smaller. When neutrinos come mainly 
from charged pion decays, « 2, we observe basically 
two regions: below sin^ i?23 — 0.52, the flavor ratio at the 
Earth lies around 0.5, whereas for larger angles it grows 
until values around 0.6. This is a consequence of the flat 
behavior of A at larger than one for z?23 in the first 
octant, see the upper panel of Fig.|5| This difference can 
even increase at the bump, where i?^ can reach values 
larger than 0.7 for -§23 in the second octant whereas it 
remains smaller than 0.6 for angles in the first octant. 
Finally at higher energies (charm neutrinos) it becomes 
very difficult to disentangle the different values: i?° ~ f 
and Rp becomes roughly 0.45 for all 7923 around 7r/4, see 
upper panel of Fig. O 




12 3 4 




12 3 4 

R" 

FIG. 9; Contours of the coefficients A{R° ,■023) (top) and 
B{R°,-Q2i) (bottom) in terms of i?° and sin^i?23. 

Let us now briefiy discuss the dependence of i?^ on diz 
and (Sep. According to Eq. Q, this variation depends 
on i?" and 1923 through the function B{R'^ ,^23)- The 
bottom panel of Fig shows that the maximal variation 
takes place when 7^23 is in the first octant and i?° deviates 
significantly from one. In all cases, though, the sensitiv- 
ity does not exceed 10%. In the bottom panel of Fig. 1101 
we show the variation of R for values of 1913 between 
and 10^^, assuming Sep = and sin^ i923 = 0.35 (green) 
and 0.65 (yellow). While at low energies (i?° ~ 0) both 
cases give similar results, at higher R^ the dependence 
on 79i3 is almost negligible. We show also the variation 
of i?p as function of 5cp for the same values of sin^ 7923 
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and sin z?i3 = 10 ^. In this case, the variation is larger 
due to the two different signs of cos^cp- 
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fore among the most promising candidates for powerful 
neutrino sources. We have calculated the yield of high 
energy neutrinos from sources with photons as target ma- 
terial, paying especially attention to opaque sources. We 
have found that the neutrino spectra from meson and 
muon decays are strongly modified with respect to trans- 
parent sources as soon as multiple scattering of nucleons 
becomes important. The main consequence is a strong 
suppression of the neutrino flux from pion and kaon de- 
cays at energies above a critical energy E^, defined as the 
energy at which the scattering length of mesons equals 
the decay length. Above this energy both pions and 
kaons scatter before decay and thus the main contribu- 
tion to the neutrino flux is supplied by the decay of charm 
mesons. As a result, the parameter range where opaque 
sources can produce ultra-high energy neutrinos is rather 
restricted. An opaque source with Ecr > 10^^ eV requires 
basically a large extension with low densities. 

Both neutrino telescopes and extensive air shower ex- 
periments can not only detect high energy neutrinos but 
also they have some flavor discrimination possibilities. 
Since the determination of the neutrino flavor ratio can 
shed light on the properties of both sources and neutri- 
nos, we have calculated also the expected flavor ratio for 
opaque sources containing photons as scattering targets. 
The main characteristic is a strong energy-dependence of 
the ratio i?° of and Ve fluxes at the sources. A generic 
prediction for neutrinos from hadron-photon scattering 
is the flavor ratio i?° below the threshold, i.e. a flux 
dominated by 9^ at the source. In the case of opaque 
sources the flavor ratio presents a bump close to the en- 
ergy where the decay length of charged pions and kaons 
equals their interaction length on target photons. The 
neutrino spectra at the source are modulated by oscilla- 
tion. Therefore the expected flavor ratio i?^ at the Earth 
will be different from the original one i?'^: it signiflcantly 
depends on the neutrino mixing parameters, and espe- 
cially on '023- Therefore the observation of a strong en- 
ergy dependence of ii^ will not only be a hint on the 
kind of source but also allow obtaining complementary 
information on the neutrino mixing angles. 



FIG. 10: (Color online) Top: Flavor ratio at the Earth i?^ 
for different values of sin^ ■!?23 assuming sin'^ i9i2 — it/6 and 
^13 = for a source at T = 10^ K and interaction depth 
r = 40. Bottom: Dependence of on i9i3 for Sep = in 
yellow (very light gray) and green (light green), and on Sep, 

assuming sin^ -ffn = 10"^ in red (gray) and blue (dark gray). Acknowledgments 
for two extreme values of sin^ i923. 



IV. SUMMARY 

The high-energy cosmic ray flux from opaque or hidden 
neutrino sources is suppressed and these source are there- 
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APPENDIX: TWO TWO-COLUMN FORMULAE 

The explicit expression for A and B are 

. ^ 12 + 39i?P - 12(fl° - 1) cos(2t?23) + 13fl° cos(4t?23) ..^ 
^ ' 52 + 25i?o + 12(i?0-l)cos(2i?23)-13i?Ocos(4i?23) 

„.^0 a ^ ^ -512x/3(l + i?")(-l + + cos(2t^23)) sin(2^23) , . 

^ ' [52 + 25i?" + 12(i?0-l)cos(2?923)-13i?0cos(4i?23F' ^ ' 
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